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1  Introduction 


This  is  the  second  in  a  series  of  three  papers  dealing  with  the  time  rates  of 
generalized  strain  tensors.  The  main  results  of  Part  I  (Scheidler  [6])  may  be 
summarized  as  follows.  Let  U  =  A,!!,®!!,  and  V  =  A,Vi®v;  denote 

the  right  and  left  stretch  tensors,  respectively,  corresponding  to  a  C2  motion. 
Here  {A,}  are  the  principal  stretches,  {u,}  is  any  principal  basis  of  U,  and 
{v,}  is  the  corresponding  principal  basis  of  V.  Let  R+  denote  the  positive 
reeds.  Then  for  any  C1  function  /  :  R+  — ♦  R,  the  material  time  derivative  of 
f(U)  =  /(A,)u,  ®  u,  and  the  Jaumann  rate  of  f(V)  =  ]r^=1  /(AJv,  ®  v, 

are  given  by  the  component  formulas 

3 

^  F(A„AJ)Z)1Ju1  ®  Uj .  (1.1) 

>0=1 

53  F(At.  ©  v_, .  (1.2) 

>o=i 

where 


f(U)'  = 
f(V)°  - 


F(A..  Aj) 


F(X,.XJ)  =  A, /'(A,) 

2A, Aj 

A(  4-  A j  A,  Xj 

y  +  y  /iK]  -/(>,) 


A  |  I  A  j  A  |  A  j 


if  A,  =  \j , 

if  A,  ^  A_, , 
if  A,  ^  \J , 


(1.3) 

(1.4) 

(1-5) 


and  {D,j}  are  the  components  of  the  stretching  tensor  D  relative  to  {v,}. 
Recall  that  the  Jaumann  rate  f(V)°  is  related  to  the  material  time  derivative 
f(Vy  by 

f(V)‘  =  f(V)°  +  Wf(V)  -  f(V)W  ,  (1.6) 

where  W  denotes  the  spin  tensor.  When  /( 1)  =  0,  /'( 1)  =  1  and  /'  >  0. 
/  is  called  a  strain  measure ,  and  f(U)  and  f(V)  are  called  the  generalized 
Lagrangian  and  Eulenan  strain  tensors ,  respectively,  corresponding  to  the 
strain  measure  /. 

Because  of  their  simplicity  these  formulas  are  useful  in  theoretical  analy¬ 
ses,  as  demonstrated  for  the  formula  (1.1)  by  Hill  [3,4].  On  the  other  hand, 
to  actually  compute  f(V)°,  say,  for  a  given  motion  from  the  formulas  above 
requires  the  calculation  of  the  eigenvalues  and  eigenvectors  of  V  at  each  place 
and  time.  Henre  it  would  also  be  useful  to  have  a  simple  expression  for  the 
tensor  f(V)°  directly  in  terms  of  the  tensors  V  and  D,  or  in  terms  of  the  left 
Cauchy-Green  tensor  B  and  D.  Unfortunately,  as  we  will  show  in  Part  III. 
the  coefficients  in  such  basis-free  formulas  are  extremely  complicated  when 
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the  principal  stretches  are  distinct.  Thus  we  are  led  to  seek  simple  approx¬ 
imations  to  the  basis-free  formulas  for  f(U)’  and  f(V)°.  In  principle,  these 
approximate  formulas  could  be  derived  from  the  exact  basis-free  formulas 
obtained  in  Part  III.  However,  the  complexity  of  the  exact  basis-free  formu¬ 
las  makes  this  approach  impractical.  Instead,  we  will  derive  the  approximate 
basis-free  formulas  from  the  exact  component  formulas  above.  Our  approxi¬ 
mate  formulas  involve  an  arbitrary  parameter  A  which  can  always  be  chosen 
in  such  a  way  that  the  formulas  provide  good  estimates  when  the  shear  strains 
are  small,  regardless  of  the  volumetric  strain  and  the  strain  rates. 

In  Section  2  we  introduce  some  notation  and  state  the  assumptions  used 
in  deriving  the  approximate  formulas.  In  Section  3  we  present  the  main 
results  of  this  paper,  namely  the  approximate  basis-free  formulas  for  f(U)' 
and  f(V)°,  together  with  explicit  bounds  for  the  error  in  these  formulas. 
In  Section  4  we  list  some  approximate  formulas  involving  the  spin  O  of  the 
rotation  tensor,  as  well  as  some  approximate  formulas  for  the  tensors  f(U) 
and  f(V).  The  derivation  of  the  results  in  Sections  3  and  4  is  sketched  in 
Section  5.  In  Section  6  we  apply  the  general  results  of  the  preceding  sections 
to  the  logarithmic  strain  tensors.  We  give  a  rigorous  proof  of  an  approximate 
formula  for  (InU)'  due  to  Hill  [4],  and  we  obtain  an  improved  version  of  an 
approximate  formula  for  (In  V)°  due  to  Gurtin  and  Spear  [2], 

2  Notation  and  basic  assumptions 

The  derivation  of  the  approximate  formulas  involves  the  use  of  Taylor's  The¬ 
orem  to  express  /(A,)  and  /'(A,)  (?  =  1.2.3)  in  terms  of  /  and  its  derivatives 
evaluated  at  some  number  A.  Our  results  take  the  general  form 

r  =  rn(A)  +  £"A„,  (2.ii 

where  n  is  a  positive  integer  (usually  1  or  2).  T„(A)  is  the  basis-free  approx¬ 
imation  to  some  tensor  T  such  as  f(U)‘  or  f(V)°,  and 

s  =  max  | A,  -  A| .  (2.2) 

»=1,2,3 

We  determine  explicit  bounds  on  the  norm  of  the  term  An.  This  allows  us  to 
bound  the  error  in  approximating  T  by  T„(A).  Roughly  speaking,  F„(A)  is  a 
good  approximation  to  T  when  e  is  sufficiently  small,  i.e.,  when  all  principal 
stretches  are  sufficiently  close  to  A.  Thus  for  a  given  set  of  principal  stretches 
a  reasonable  choice  for  A  is  the  one  which  minimizes  e.  However,  in  view  of 
the  dependence  of  Tn  on  A,  the  most  useful  choice  for  A  will  generally  depend 
on  the  intended  applications  of  the  approximate  formula;  some  examples  are 
given  at  the  end  of  this  section.  Unless  stated  otherwise,  we  regard  A  as 
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arbitrary.  Note  that  since  the  principal  stretches  are  time-dependent  scalar 
fields,  in  general  A,  e  and  the  other  scalars  introduced  below  are  also  time- 
dependent  scalar  fields. 

Let 

^min  =  min{A! ,  A2,  A3} ,  Amax  =  max{A1?  A2,  A3) .  (2.3) 

Then 

e  =  max{AmflX  -  A,  A  -  Am,„}  .  (2.4) 

For  any  symmetric  tensor  A  let  5(A)  denote  the  maximum  orthogonal  shear 
component  of  A: 

s(A)  =  max(a  •  Ab) ,  (2.5) 

a.b 

where  a  and  b  range  over  all  pairs  of  orthogonal  unit  vectors.  Then  (Erick- 
sen  [1,  §46]) 

-S(U)  =  s(V)  =  (Amoi  -  Amm)/2  =  s .  (2.6) 

Thus  s  is  the  largest  shear  strain  corresponding  to  the  stretch  tensors  U  and 
V.  From  (2.4)  it  follows  that 

s  <  £  (2.7) 

for  all  choices  of  A,  and 

£<2s.  if  Amin  <  A  <  Amax .  (2.S) 

Hence,  if  £  is  small  then  the  shear  strains  are  small.  Conversely,  if  A  satisfies 
Amin  <  A  <  Amax  and  if  the  shear  strains  are  small,  then  £  is  small.  Also 
note  that 

£<||U  — AI||  =  ||V  — AI||<  y/Zs.  (2.9) 

Here  ||  •  ||  denotes  the  norm  of  a  tensor;  i.e. ,  for  any  tensor  H, 

3  1,/2 

IIHII  £  (tr(HTH)]1/2  =  'EM,)*  .  (2.10) 

>J=1 

where  {HtJ}  are  the  components  of  H  relative  to  any  orthonormal  basis.  The 
inequality  (2.9)  implies  that  £  is  small  iff  U  and  V  axe  close  to  the  dilatation 
AI.  All  the  results  in  this  paper  axe  valid  under  the  following  conditions: 

1.  £  <  A;  in  view  of  (2.4)  this  is  equivalent  to  A  >  Amor/2. 

2.  The  motion  is  C2. 

3.  The  function  /  :  R+  — ♦  R  is  C3  A 

'The  one  exception  is  equation  (4  6),  which  involves  the  first,  k  derivatives  of  /  for 
arbitrary  k.  However,  only  the  cases  k  <  3  will  be  utilized  in  the  derivation  of  other 
formulas. 


In  particular,  we  do  not  require  that  /  be  a  strain  measure,  although  this  is 
the  case  of  most  interest.  The  condition  £  <  A  is  imposed  solely  to  obtain 
simple  error  bounds.  For  some  choices  of  A  this  condition  is  automatically 
satisfied,  while  for  other  choices  of  A  it  imposes  a  constraint  on  the  principal 
stretches  (see  the  examples  below). 

The  condition  e  <  A  is  equivalent  to  the  existence  of  a  scalar  field  r  such 
that 

£  <  r  A,  0  <  r  <  1 .  (2.11) 

This,  together  with  the  definition  of  e ,  implies  that 

A,  €  /r(A)  =  [(1  —  r) A,  (1  +  r)A]  C  R+  ,  *  =  1,2,3.  (2.12) 

The  factor  1/(1  —  r)  will  appear  in  some  of  the  error  bounds;  hence  these 
bounds  are  useful  only  for  those  motions  for  which  r  is  bounded  away  from 
1.  Some  of  our  results  take  a  simpler  form  when  expressed  in  terms  of  the 
scalar 


i  =  max  — -  —  1  . 

.=  1,2,3  A 

(2.13) 

By  definition 

(2.2)  we  have 

£  —  Xi . 

(2.14) 

so  that  (2.11) 

is  equivalent  to  the  condition 

i  <  r  <  1  . 

(2.15) 

In  the  remainder  of  this  section  we  discuss  some  appropriate  choices  for  A. 
Example  I.  By  (2.4)  we  see  that  the  value  of  A  which  minimizes  e  and 

i  is 

A  =  (Amin  +  Amar)/2 ,  (2.16) 


in  which  case  e  =  (Amai  ~^m>n)/2  =  s.  Thus  £  is  small  iff  the  shear  strains  are 
small.  Note  that  the  condition  e  <  A  is  satisfied  for  all  motions.  Equivalently, 
(2.11)  holds  for 


^mfli 

Amln 


(2.17) 


in  which  case 


(2.18) 


Thus  1/(1  -  r)  is  bounded  iff  the  motion  is  such  that  Amoi/Amin  is  bounded. 

Example  II.  At  the  other  extreme,  we  may  simply  choose  A  =  1.  When 
/  is  a  strain  measure  this  choice  yields  the  simplest  expressions  for  the  co¬ 
efficients  in  the  approximate  formulas.  As  we  will  see  in  §6.  this  choice  is 
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sometimes  useful  for  comparison  with  other  approximate  formulas  in  the  lit¬ 
erature.  Note  that  the  bounds  Amin  <  A  <  Aroox  need  not  hold  in  this  case, 
and  that  the  condition  e  <  A  reduces  to  e  <  1,  which  would  usually  be 
imposed  anyway.  Also  note  that  e  is  small  iff  all  principal  stretches  are  suffi¬ 
ciently  close  to  1  iff  all  principal  strains  are  sufficiently  close  to  zero.  Hence 
for  this  choice  of  A  not  only  the  shear  strains  but  also  the  volumetric  strain 
must  small  be  to  obtain  good  approximations. 

Example  III.  A  natural  choice  for  A  is  the  geometric  mean  of  the  prin¬ 
cipal  stretches: 

A  =  (A1A2A3)1/3.  (2.19) 

Then 


A  =  (det  U)1/3  =  (det  V)’/3  =  (det  F)1/3  = 


(2.20) 


where  F  denotes  the  deformation  gradient,  p  and  v  denote  the  mass  density 
and  the  specific  volume,  respectively,  in  the  deformed  configuration,  and  p0 
and  v0  denote  the  corresponding  quantities  in  the  reference  configuration. 
Thus  A  may  also  be  interpreted  as  the  volumetric  stretch.  Unlike  (2.16).  this 
choice  of  A  does  not  require  knowledge  of  the  principal  stretches  since  A  can 
be  determined  directly  from  (2.20).  In  particular,  A  =  1  for  any  isochoric 
motion.  Since  (2.19)  implies  Amin  <  A  <  \max,  it  follows  from  (2.7)  and 
(2.S)  that  s  <  £  <  2s.  Thus  £  is  small  iff  the  shear  strains  are  small.  Note 
that  the  condition  £  <  A  does  not  hold  in  general:  e.g.,  if  A,  =  A2/2  and 
A3  =  2A2.  then  e  —  A  =  A2.  However,  since  a  small  £  is  required  for  good 
approximations,  for  many  applications  the  constraint  £  <  A,  or  equivalently, 
the  condition  (2.11),  is  fairly  mild.  For  example,  in  many  solids  elastic  shear 
strains  cannot  exceed  19c,  and  density  increases  by  a  factor  greater  than  2 
are  extremely  difficult  to  achieve.  Hence,  for  elastic  deformations  of  these 
solids  we  may  take  s  <  0.01  and  p/p0  <  2,  in  which  case  £  <  2s  <  0.02  and. 
by  (2.20).  A  >  0.79.  so  that  (2.11)  is  satisfied  for  r  =  0.03. 

For  A  given  by  (2.19)  it  is  useful  to  introduce  the  right  and  left  distortional 
stretch  tensors  U  and  V,  respectively: 

v  =  \v,  V  =  i  V.  (2.21) 

A  A 


By  (2.20)  it  follows  that  U  and  V  are  unaffected  by  the  dilatational  part 
of  the  deformation  in  the  sense  that  the  deformation  gradients  F  and  aF 
have  the  same  distortional  stretch  tensors  for  any  a  >  0;  in  particular, 
det  U  =  det  V  =  1.  The  eigenvalues  { A,  / A }  of  U  and  V,  and  the  eigen¬ 
values  {A, /A  —  1}  of  U  —  I  and  V  —  I.  are  the  principal  distortional  stretches 
and  strains,  respectively.  These  scalars,  and  hence  £,  are  also  unaffected  by 
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the  dilatational  part  of  the  deformation.  The  condition  e  <  A,  or  equiv¬ 
alently.  the  condition  (2.15),  is  simply  the  requirement  that  the  principal 
distortional  strains  have  magnitude  less  than  1. 

Example  IV.  There  are  other  choices  for  A  which  retain  the  interpreta¬ 
tion  of  a  mean  stretch.  Examples  are 

X  =  \(  Ax  +  A2  +  A3)  =  ~trU  =  itr  V  (2.22) 

and 

»  =  +  +  *ll  'n  =  ^  (tr  C)I,J  =  2=  (tr  B>l/S,  <’23) 

where  C  =  FTF  =  U2  and  B  =  FFT  =  V2  are  the  right  and  left  Cauchy- 
Green  tensors. 

3  Main  results 

Let  R  denote  the  rotation  tensor  and  let  D  denote  the  rotated  stretching 
tensor:  D  =  RTDR  Let 

a(A)  =  l-[f[  A)  +  A /"(A)]  =  ~[A/'(A)1 .  (3.1) 

In  particular.  a(A)  =  (1  +  /"(l))/2  when  A  =  1  and  /  is  a  strain  measure. 
We  will  show  that 

f(U)-  =  A /'(A)  D  +  c 

=  -A2/"(A)  D  +  q(A)(DU  +  UD)  +  e2$2  -  (3.2) 

and 

f(V)°  =  A/'(A)D  +  £*3 

=  -A2/"{A)D+a(A)(DV+VD)  +  e2$4.  (3.3) 

Here  and  below,  (k  =  1,2.. ..,8)  is  a  symmetiic  tensor  satisfying 

lhM<**(A)||D||  (3.4) 

for  some  continuous  function  <pk.  When  the  terms  of  order  e  or  e2  are  ne¬ 
glected,  equations  (3.1)— (3.3)  yield  approximate  basis-free  formulas  for  f(U)‘ 
and  f(V)°. 

For  any  tensor  field  H,  we  have 
H  f  ( V)  —  f(V)H  = 

=  /'(A)(HV  -  VH)  +  •  (3.5) 


6 


Here  and  below,  k  (k  =  1,2,3)  is  a  tensor  satisfying 

11**11  <V>*(A)||H||  (3.6) 

for  some  continuous  function  xl'k.  Then  (3.5)  with  H  =  W,  together  with 

(3.3)  and  (1.6),  yield  approximate  basis-free  formulas  for  f(V)‘. 

Let  O  denote  a  tensor- valued  function  of  a  nonnegative  real  variable  such 
that  ||(9(/i)||  <  Mh  for  some  constant  M  and  all  sufficiently  small  h  >  0. 
Then  (3.3),  and  (3.4)  imply  that  f(V)°  =  A/'(A)D  +  O(e).  Similarly,  we 
could  replace  the  term  £2$4  in  (3.3)2  by  0(e2).  This  is  typically  the  way 
that  approximate  formulas  are  stated  in  the  literature;  cf.  the  examples  in 
§6.  While  such  results  are  of  theoretical  interest,  they  axe  of  little  practical 
use.  If  we  know  onlv  that  f(V)°  =  r„  +  0(£n),  then  we  may  conclude  that 
l|f(V)°  -  rj|/e«-*  -  0  as  £  — ►  0,  but  for  a  given  value  of  e  we  cannot 
determine  the  error  in  approximating  f(V)°  by  Tn.  To  do  this  we  need 
explicit  bounds  for  the  remainder  terms.  These  bounds  can  be  obtained  from 

(3.4)  and  (3.6)  provided  that  we  have  explicit  for  formulas  for  the  functions 
spk  and  vv  Such  formulas  are  given  below. 

The  stretch  tensors  U  and  V  are  irrational  functions  of  the  deformation 
gradient  F.  For  computational  purposes  it  would  be  useful  to  have  formulas 
analogous  to  (3.2)2.  ( 3.3)2  and  (3.5)2.  but  involving  the  Cauchy-Green  tensors 
C  and  B  instead.  Let 


$(A)  =  ~[/'(A)  -  A/"(A)] . 


‘  A 


+  /"(A) 


o(A) 

2A 


(3.7) 


(3.S) 


In  particular,  3{ A)  =  (1  -  /"( 1 ) )/2  and  7(A)  =  (1  4-  /"(l))/4  when  A  =  1 
and  /  is  a  strain  measure.  The  analogs  of  (3.2)2,  (3.3)2  and  (3.5)2  are 


f(U)'  =  /?(A)D  +  7(A)(DC  +  CD)  +  £2*5,  (3.9) 

f(V)°  =  /?(A)D  +  7(A)(DB  +  BD)  +  £2$6,  (3.10) 

and 

Hf(V)-f(V)H  =  ^(HB-BH)  +  e2$3  (3.11) 

Z  A 

for  any  tensor  field  H. 

Now  suppose  that  /'(A)  ^  0,  as  is  the  case  when  /  is  a  strain  measure. 
Let 


/'( A)  = 


1 

2 


1  + 


A /"(A) 


/'(A)  j 


t'(A)  =  A /'(A)  -  2/i(A)/(A) . 


(3.12) 


In  particular.  p( A)  =  (1  +  /"(l))/2  and  i/(\)  =  1  when  A  =  1  and  /  is  a 
strain  measure.  We  will  show  that 

f(U)'  =  j/(A)D  +  /i(A)[Df(U)  +  f(U)D)  +  e2$7  (3.13) 

and 

f(V)‘  =  i/(A)D  +  [//(A)D  +  W]f(V)  +  f(V)[/z(A)D-W] 

+£2*8-  (3-14) 

When  the  terms  of  order  e2  cure  neglected,  (3.13)  and  (3.14)  yield  approxi¬ 
mate  linear  ordinary  differential  equations  for  f(U)  and  f(V).  For  a  given 
material  point  X,  integration  of  (3.14)  (with  neglected)  requires  only 
the  time  histories  of  A  and  the  stretching  and  spin  tensors  corresponding  to 
X.  together  with  the  value  of  f(V(X,f0))  at  some  initial  time  t0.  Note  that  if 
A  is  taken  to  be  the  volumetric  stretch  (cf.  2.20),  then  A  satisfies  the  ordinary 
differential  equation 

A  =  -(tr  D)A  .  (3.15) 

J 

To  obtain  bounds  for  the  error  when  the  remainder  terms  are  neglected 
in  the  formulas  above,  we  need  explicit  expressions  for  the  functions  <pk  and 
Vic  in  (3.4)  and  (3.6).  Recall  the  definition  (2.12)  of  the  interval  /r(A).  Let 
f(k)  denote  the  kth  derivative  of  /,  and  let 

/*(A)=  max  |/(fc)(x )| .  (3.16) 

x€/r(A) 

Then  we  may  take 


v,(A) 
^’3  ( A ) 

*i(A) 


2/i  (A),  V>2(A)  =  /2(A), 

I/'(A)|  +  2A/3(A) 
V»3(A)  = - : - 


<p2{\)  =  94(A)  = 


l/'(A)| 


+  l/"(A)|  +  A/3(A) 


V’s(A)  =  vp6(A)  =  3v?2(A)/2, 

i^7(A)  =  ^(A)  =  <r,2(A)  +  2|/r(A)|/2(A) . 


(3.17) 
(3. IS) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 


If  f"  has  constant  sign,  then  /2(A)  may  be  replaced  by  f2( A)/2  in  (3.17)2, 
(3.18)  and  (3.19).  If  f"  has  constant  sign,  then  /3(A)  may  be  replaced  by 
/.a ( A ) / 2  in  (3.20).  In  view  of  (2.14),  we  may  replace  e  by  i  in  the  remainder 
terms,  provided  that  the  right-hand  side  of  (3.17)]  and  (3.19)  is  multiplied  by 


S 


A,  and  the  right-hand  side  of  (3.17)2.  (3.18)  and  (3.20)-(3.22)  is  multiplied 
by  A2. 

For  the  special  case  f(x)  =  x  we  have  f(V)  =  V,  and  the  results  in  this 
section  reduce  to 

V°  =  AD +£$3 

=  ^(DV  +  VD)  +  £2$4 

=  ^D  +  ^DB  +  BDJ  +  e**.,  (3.23) 

and 

V  =  V°  +  WV-VW 

=  V  +  (3.24) 

where 

ll*j||<^l|W||, 

ll*3l|.  *ll*«l|.  ^Il*6ll  <7^-. 

3  1  —  r 

4  Some  additional  formulas 

Let  Cl  denote  the  spin  of  the  rotation  tensor:  Cl  =  RRT  =  —  ClT .  In  addition 
to  the  Jaumann  rate  of  f(V).  the  corotational  rate  f(V)"  is  also  of  interest. 
Recall  (Part  I,  §5)  that  f(V)"  is  defined  in  terms  of  O  and  the  material  time 
derivative  of  f(V)  by 

f(V)'  =  f(V)‘  +  flf(V)  -  f(V)fi.  (4.1) 

Approximate  formulas  for  the  term  fif(V)  —  f(V)fi,  as  well  as  bounds  for 
the  error  in  these  formulas,  are  obtained  by  setting  H  =  Cl  in  (3.5),  (3.6)  and 
(3.11).  Moreover,  we  will  show  that  f(V)°  may  be  replaced  by  f(V)’  in  (3.3) 
and  (3.10).  Of  course,  the  tensors  $4  and  $6  in  the  remainder  terms  will 
generally  be  different,  but  they  satisfy  the  same  bounds  as  before.  It  follows 
that  f(V)*  and  f(V)°  agree  to  within  terms  of  order  £2.  In  fact,  we  have 

f(V)*  =  f(V)°  +  .  (4.2) 

Here  and  below.  $k  (k  =  1,2, 3, 4)  is  a  symmetric  tensor  satisfying 

ll*dl  <  ft(A)||D||  (4.3) 


(3.25) 

(3.26) 
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for  some  continuous  function  (pk . 

Since  equations  (3.2),  (3.9)  and  (3.13)  for  f(U)‘  involve  the  rotation  tensor 
R  through  the  rotated  stretching  tensor  D,  and  since  R  =  fiR,  approximate 
formulas  for  $7  in  terms  of  W  are  also  of  interest.  We  have 


ft  —  W  +  £  4»2 

=  W+^(DV-VD)-he2$3 

=  W  +  ±(DB-BD)  +  e2*4. 

For  the  functions  ^>k  in  the  above  formulas  we  may  take 

.  2A/j  (A)  .  4.  1 

V'1-(i_r)3'  ^2-V>3-5S?4-1_r 


(4.4) 


(4.5) 


Finally,  we  consider  some  approximate  formulas  for  f(V)  and  f(U).  If  / 
is  Ck  (k  >  1)  then 


f(V)  =  E— “(V-Air  +  e^ 

n : 


n= 0 


(4.6) 


where  Afc  is  a  sj-mmetric  tensor  satisfying 


l|Ajt||  <  —  A(A) . 

Recall  that  fk  is  defined  by  (3.16).  Let 


(4.7) 


VW  =  /"(A)  - 


/'(A) 


(4.S) 


Then  we  also  have 


/'(A), 


f(V)  =  /(A)  I+^fB-A'IJ  +  e’A,, 

=  /(-')I+^(B-A!I)  +  ^(B-AU)’+eA3,  (4.9) 


where  A2  and  A3  are  symmetric  tensors  satisfying 


|A2||<  jil^!  +  /2(A),  ||A3||<1 


PMUI 


+  /a(A) 


(4.10) 


To  obtain  formulas  for  f(U),  simply  replace  V  by  U  and  B  by  C  in  (4.6)  and 
(4.9).  Note  that  77 ( A )  =  /"( 1)  —  1  when  A  =  1  and  /  is  a  strain  measure. 
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5  Derivation  of  the  results 


Equations  (3.1)— (3.3),  (3.4)  for  k  <  4,  (3.19)  and  (3.20)  follow  from  (1.1)— 
(1.5)  and  Taylor’s  Theorem  for  /  and  /'  with  Lagrange’s  formula  for  the 
remainder.  In  deriving  (3.2)  we  have  used  the  fact  that  {DtJ}  are  also  the 
components  of  D  relative  to  {u,}.  Id  deriving  (3.20)  we  have  also  used  the 
following  result  due  to  Hummel  and  Seebeck  [5]: 


f(y)-f(x)  f'(y)  +  f'(x)  f"'(6) 


for  some  0  between  x  and  y. 

We  illustrate  the  method  by  deriving  (3.3)2.  We  must  show  that 

F(  A,.  A,)  -  —  A2/"(  A)  +  i[/'( A)  +  A/"(A)](  A,  +  A,)  +  6ijS*  .  (5.2) 

where 

l^.jl  <  <r'4 ( A )  (5.3) 

for  v4(A)  satisfying  (3.20)2.  Then  substitution  of  (5.2)  into  (1.2)  yields  ( 3.3 )2 
with  q(A)  given  by  (3.1)  and 


3 

$4  =  H  t.jD.jV,  ©  V, .  (5.4) 

<0=1 

Then  (3.4)  with  k  =  4  follows  from  (5.4)  and  (5.3).  Thus  it  remains  to 
establish  (5.2)  and  (5.3).  We  will  use  the  following  inequalities,  which  are 
consequences  of  (2.2),  (2.11)  and  (2.12): 


I  A,  -  Xj\  <  2s  ,  (5.5) 

2(1  -  r) A  <  A,-  +  A,  <  2A  +  2e  <  2(1  +  r) A  ,  (5.6) 

e2  <  r\e  <  r2 A2  .  (5.7) 

Suppose  that  A,  ^  \y  Then  F(A,,  \})  is  given  by  (1.5).  By  (5.5)  and  (5.6)j, 


A2  +  A^ 

_j _ i 

A,  +  A; 


P>} 


A,  +  Aj 

2 


+  P.j 


A  +  -(A,  +  Xj  -  2A)  +  p,} , 

(5.8) 

(A,  —  Aj)2  ___  e2 

2(A,  +  X})  ~  (1  -r)A  ' 

(5.9) 
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By  applying  the  Tajdor- Lagrange  formula  to  /'  we  obtain 

f'(K)  =  /'(A)  +  /"(A)(A.  -  A)  +  i/"'(2,)( A,  -  A)!  (5.10) 

for  some  number  z,  between  A,  and  A,  By  setting  y  —  A,  and  x  =  in  (5.1), 
and  using  (5.10)  and  an  analogous  formula  for  /'( we  obtain 

=  /'(A)  +  +  A,  -  2A)  +  q,}  ,  (5.11) 

%  =  -  A)2  +  A,  -  A)2  -  ^rm A,  -  A,)2.  (5.12) 

Since  Zj  and  #  all  lie  in  the  interval  /(A),  by  (5.12),  (5.5),  (3.16)  and  (2.2) 
we  have 

l?J</3(A)£2.  (5.13) 

From  (1.5),  (5.8)  and  (5.11)  we  obtain  (5.2)  with 

v2  -  Pan A)+  5(A,  +  AJ-2A)  +  ptJ  .I/"(A)(A,  +  AJ-2A) 

+  ~(A,  d-AjJ  +  p.j  qv  .  (5.14) 

Then  (5.3).  with  ip4(A)  given  by  (3.20)2,  follows  from  (5.14),  (5.13).  (5.9). 
(5.7)  and  (5.6).  A  similar  analysis  shows  that  (5.2)  and  (5  3)  also  hold  when 
A,  =  Xj.  This  completes  the  proof  of  (3.3)2. 

Equations  (3.5),  (3.6)  for  k  <  2,  and  (3.17)  follow  from 

H  f(V)  -  f(V)H  =  ^[/(^-/(AJlf^v.©^,  (5.15) 

*0=1 

and  the  Taylor- Lagrange  formula  for  /. 

Equations  (3.7)— (3.1 1 ),  (3.4)  for  k  =  5  and  6,  (3.6)  for  k  =  3,  (3.18) 
and  (3.21)  follow  from  the  corresponding  results  in  terms  of  U  and  V,  the 
identity 

v  =  5i  +  Sb-JT(V-ai)j’  (5-16) 

and  an  analogous  identity  involving  U  and  C.  Equations  (3.9)  and  (3.10)  can 
also  be  derived  without  recourse  to  (3.2)  and  (3.3).  For  example,  (3.10)  can 
be  obtained  directly  from  the  component  formula  I-(5.20)2  with  h(x 2)  =  f(x). 

2An  equation  number  prefixed  by  “I"  denotes  an  equation  in  Part  I. 
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Equations  (4.6)  and  (4.7)  follow  from  f(V)  =  £,3_i  /( A,)v,  0  v,  and  the 
Taylor- Lagrange  formula  for  /.  Equations  (4.8) — (4.10)  follow  from  (4.6) 
and  (4.7)  (with  k  =  2,3),  the  identity  (5.16),  and  (2.9). 

By  solving  (4.6)  with  k  =  2  for  V,  substituting  the  result  into  (3.3)2,  and 
using  (3.4),  (4.7)  and  (1.6),  we  obtain  (3.14),  (3.12),  (3.4)  for  k  =  8,  and 
(3.22)2.  Equations  (3.13)  and  (3.22)!  are  derived  similarly. 

To  see  that  f(V)°  may  be  replaced  by  f(V)*  in  (3.3)  and  (3.10),  use 
f(V)'  =  Rf(U)  RT  (cf.  I-(5.9)),  (3.2),  (3.9)  and  ||R$*RT||  =  j|$J.  Equa¬ 
tions  (4.2),  (4.3)  for  k  —  1,  and  (4.5) j  follow  from  the  component  formula 
I-(5.16).  Equations  (4.4),  (4.5)  and  (4.3)  for  k  >  2  follow  from  the  component 
formula  I- ( 5. 14)  and  the  identity  (5.16). 


6  The  logarithmic  strain  tensors 

Approximate  basis-free  formulas  for  (InU)’.  (In  V)°  and  (In  V)’  are  obtained 
from  (3.1).  (3.2)2.  (3.3)2,  (3.4).  (3.20).  the  comments  following  (3.22).  and 
the  comments  preceding  (4.2).  The  results  are 

(Ini!)'  =  D  +  f  20,  (6.1) 


and 

(In  V)°  =  D  +  £  2©2 .  (InV)*  =  D  +  e2©3. 
where  the  symmetric  tensors  satisfy 

3 


l©*ll  < 


(1  -r)4 


i  D  | 


k  =  1.2.3. 


(6.2) 


(6.3) 


From  (3.5),  (3.6),  (3.11).  (3.17),  (3.18)  and  the  comments  following  (3.22). 
we  obtain  the  following  results.  For  any  tensor  field  H. 

H(lnV)  -  (InV)H  = 

=  Y(HV~VH)  +  f~2*2 
A 


where  the  tensors  ’d'*  satisfy 


for 


Vi  = 


1  —  r 


-(HB-BH)  +  f2$3, 

(6.4) 

ll<  F  =  1,2,3, 

(6.5) 

1 

(6.6) 
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Since 


(In  V)'  =  (In  V)°  4-  W(ln  V)  —  (In  V)W 

=  (InV)*  4-  ft(lnV)  -  (lnV)fi,  (6.7) 

approximate  basis-free  formulas  for  (InV)'  follow  from  (6.2)  and  (6.4)  with 

H  =  W  or  H  =  fi. 

The  approximate  formula  (6.1)  is  essentially  due  to  Hill  [3,4],  although 
his  derivation  is  generally  valid  only  for  the  case  of  distinct  principal  stretches. 
Hill  noted  that  for  /  =  In,  the  coefficient  F(A,.AJ)in(l.l)  satisfies  F(X,,  X})  = 

1  if  i  =  j  and  F( A,,  A^)  =  1  -  (A,/A_,  -  l)2/6  H - if  i  ±  j\  cf.  I-(6.6)  through 

I-(6.9).  He  concluded  that  (In  U)‘ =  D  +  0(||f(U)||2)  for  any  strain  measure 
/,  and  that  the  remainder  term  is  unaffected  by  the  dilatational  part  of  the 
deformation.  Hill  did  not  derive  a  bound  for  the  remainder  term. 

Gurtin  and  Spear  [2]  proved  that 

F  -  1  =  0(A)  k  F  =  O(h)  =>  (In  V)°  =  D  +  0(h3):  (6.S) 


they  did  not  derive  a  bound  for  the  remainder  term.  Their  proof  is  valid  for 
any  C2  motion  satisfying  the  conditions  on  the  left-hand  side  of  (6.S).  We 
can  easily  recover  their  result  from  ours.  Indeed,  by  setting  A  =  1  and  using 
(6.2),  (6.3),  (2.9)  and  (2.14).  we  find  that 


V  -  1  =  0(A)  fcD  =  0(A)  =► 


(In  V)°  =  D  +  0(A3) 
(In  V)*  =  D  +  0(h3). 


(6.9) 


Then  Gurtin  and  Spear's  result  follows  since,  as  shown  in  the  course  of  their 
proof,  the  conditions  on  the  left-hand  side  of  (6.8)  imply  that  V  —  I  =  0(h), 
D  =  0(h)  and  W  =  O(h).  Note  that  our  result  (6.9)  places  no  restrictions 
on  the  rotation  tensor  R  or  the  spin  tensor  W. 

Now  assume  that  A  is  defined  by  (2.16),  (2.19),  (2.22)  or  (2.23);  then  given 
any  x  >  0  there  is  some  motion  for  which  A  =  x.  Equations  (6.2)j  and  (6.3) 
imply  that  (In  V)°  =  D  +  0(i2).  We  claim  that  this  condition  characterizes 
the  logarithmic  strain  measure  in  the  sense  that  if  /  is  a  strain  measure  and 
f(V)°  =  D  4-  0(e2)  for  any  motion,  then  /  =  In.  For  by  (2.14),  (3.1),  (3.3)2 
and  (3.4),  we  see  that  this  condition  can  be  satisfied  iff  — A2/"( A)  =  1  and 
/'(A)  +  A/"(A)  =  0  for  all  A  >  0.  Since  /  is  a  strain  measure  it  follows 
that  /  =  In.  There  are  a  number  of  weaker  conditions  wffiich  also  suffice  to 
characterize  the  logarithmic  strain  measure.  By  using  (2.9)  and  arguing  as 
above,  we  see  that  if  a  strain  measure  /  satisfies  any  one  of  the  following 
conditions  then  /  =  In: 


f(V)°  =  D  +  0(||V-AI||), 

f(V)°  =  a(A)  D  4- 0(|[V  —  AI||2)  for  some  function  a  , 

f(V)°  =  D  +  6(A)(DV  +  VD)  +  0(||V-  AI|[2)  for  some  function  b. 
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In  view  of  (4.2),  we  arrive  at  the  same  conclusion  if  f(V)°  is  replaced  by 
f(V)’  in  the  above. 
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